Imprinted SARS-CoV-2 humoral immunity induces convergent Omicron RBD evolution.

Continuous evolution of Omicron has led to a rapid and simultaneous emergence of numerous variants that display growth advantages over BA.5 (ref. 1). Despite their divergent evolutionary courses, mutations on their receptor-binding domain (RBD) converge on several hotspots. The driving force and destination of such sudden convergent evolution and its effect on humoral immunity remain unclear. Here we demonstrate that these convergent mutations can cause evasion of neutralizing antibody drugs and convalescent plasma, including those from BA.5 breakthrough infection, while maintaining sufficient ACE2-binding capability. BQ.1.1.10 (BQ.1.1 + Y144del), BA.4.6.3, XBB and CH.1.1 are the most antibody-evasive strains tested. To delineate the origin of the convergent evolution, we determined the escape mutation profiles and neutralization activity of monoclonal antibodies isolated from individuals who had BA.2 and BA.5 breakthrough infections2,3. Owing to humoral immune imprinting, BA.2 and especially BA.5 breakthrough infection reduced the diversity of the neutralizing antibody binding sites and increased proportions of non-neutralizing antibody clones, which, in turn, focused humoral immune pressure and promoted convergent evolution in the RBD. Moreover, we show that the convergent RBD mutations could be accurately inferred by deep mutational scanning profiles4,5, and the evolution trends of BA.2.75 and BA.5 subvariants could be well foreseen through constructed convergent pseudovirus mutants. These results suggest that current herd immunity and BA.5 vaccine boosters may not efficiently prevent the infection of Omicron convergent variants.

SARS-CoV-2 Delta Variant: Interplay between Individual Mutations and Their Allosteric Synergy.

Since its first appearance in April 2021, B.1.617.2, also termed variant Delta, catalyzed one major worldwide wave dominating the second year of coronavirus disease 2019 (COVID-19) pandemic. Despite its quick disappearance worldwide, the strong virulence caused by a few point mutations remains an unsolved problem largely. Along with the other two sublineages, the Delta variant harbors an accumulation of Spike protein mutations, including the previously identified L452R, E484Q, and the newly emerged T478K on its receptor binding domain (RBD). We used molecular dynamics (MD) simulations, in combination with free energy perturbation (FEP) calculations, to examine the effects of two combinative mutation sets, L452R + E484Q and L452R + T478K. Our dynamic trajectories reveal an enhancement in binding affinity between mutated RBD and the common receptor protein angiotensin converting enzyme 2 (ACE2) through a net increase in the buried molecular surface area of the binary complex. This enhanced binding, mediated through Gln493, sets the same stage for all three sublineages due to the presence of L452R mutation. The other mutation component, E484Q or T478K, was found to impact the RBD-ACE2 binding and help the variant to evade several monoclonal antibodies (mAbs) in a distinct manner. Especially for L452R + T478K, synergies between mutations are mediated through a complex residual and water interaction network and further enhance its binding to ACE2. Taking together, this study demonstrates that new variants of SARS-CoV-2 accomplish both "attack" (infection) and "defense" (antibody neutralization escape) with the same "polished sword" (mutated Spike RBD).

Risk of rapid evolutionary escape from biomedical interventions targeting SARS-CoV-2 spike protein.

The spike protein receptor-binding domain (RBD) of SARS-CoV-2 is the molecular target for many vaccines and antibody-based prophylactics aimed at bringing COVID-19 under control. Such a narrow molecular focus raises the specter of viral immune evasion as a potential failure mode for these biomedical interventions. With the emergence of new strains of SARS-CoV-2 with altered transmissibility and immune evasion potential, a critical question is this: how easily can the virus escape neutralizing antibodies (nAbs) targeting the spike RBD? To answer this question, we combined an analysis of the RBD structure-function with an evolutionary modeling framework. Our structure-function analysis revealed that epitopes for RBD-targeting nAbs overlap one another substantially and can be evaded by escape mutants with ACE2 affinities comparable to the wild type, that are observed in sequence surveillance data and infect cells in vitro. This suggests that the fitness cost of nAb-evading mutations is low. We then used evolutionary modeling to predict the frequency of immune escape before and after the widespread presence of nAbs due to vaccines, passive immunization or natural immunity. Our modeling suggests that SARS-CoV-2 mutants with one or two mildly deleterious mutations are expected to exist in high numbers due to neutral genetic variation, and consequently resistance to vaccines or other prophylactics that rely on one or two antibodies for protection can develop quickly -and repeatedly- under positive selection. Predicted resistance timelines are comparable to those of the decay kinetics of nAbs raised against vaccinal or natural antigens, raising a second potential mechanism for loss of immunity in the population. Strategies for viral elimination should therefore be diversified across molecular targets and therapeutic modalities.

Cryo-EM structure of a SARS-CoV-2 omicron spike protein ectodomain.

The omicron variant of SARS-CoV-2 has been spreading rapidly across the globe. The virus-surface spike protein plays a critical role in the cell entry and immune evasion of SARS-CoV-2. Here we determined the 3.0 Å cryo-EM structure of the omicron spike protein ectodomain. In contrast to the original strain of SARS-CoV-2 where the receptor-binding domain (RBD) of the spike protein takes a mixture of open ("standing up") and closed ("lying down") conformations, the omicron spike molecules are predominantly in the open conformation, with one upright RBD ready for receptor binding. The open conformation of the omicron spike is stabilized by enhanced inter-domain and inter-subunit packing, which involves new mutations in the omicron strain. Moreover, the omicron spike has undergone extensive mutations in RBD regions where known neutralizing antibodies target, allowing the omicron variant to escape immune surveillance aimed at the original viral strain. The stable open conformation of the omicron spike sheds light on the cell entry and immune evasion mechanisms of the omicron variant.

Structural basis for SARS-CoV-2 Delta variant recognition of ACE2 receptor and broadly neutralizing antibodies.

The SARS-CoV-2 Delta variant is currently the dominant circulating strain in the world. Uncovering the structural basis of the enhanced transmission and altered immune sensitivity of Delta is particularly important. Here we present cryo-EM structures revealing two conformational states of Delta spike and S/ACE2 complex in four states. Our cryo-EM analysis suggests that RBD destabilizations lead to population shift towards the more RBD-up and S1 destabilized fusion-prone state, beneficial for engagement with ACE2 and shedding of S1. Noteworthy, we find the Delta T478K substitution plays a vital role in stabilizing and reshaping the RBM loop473-490, enhancing interaction with ACE2. Collectively, increased propensity for more RBD-up states and the affinity-enhancing T478K substitution together contribute to increased ACE2 binding, providing structural basis of rapid spread of Delta. Moreover, we identify a previously generated MAb 8D3 as a cross-variant broadly neutralizing antibody and reveal that 8D3 binding induces a large K478 side-chain orientation change, suggesting 8D3 may use an "induced-fit" mechanism to tolerate Delta T478K mutation. We also find that all five RBD-targeting MAbs tested remain effective on Delta, suggesting that Delta well preserves the neutralizing antigenic landscape in RBD. Our findings shed new lights on the pathogenicity and antibody neutralization of Delta.

Structural basis of SARS-CoV-2 Omicron immune evasion and receptor engagement.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern evades antibody-mediated immunity that comes from vaccination or infection with earlier variants due to accumulation of numerous spike mutations. To understand the Omicron antigenic shift, we determined cryo-electron microscopy and x-ray crystal structures of the spike protein and the receptor-binding domain bound to the broadly neutralizing sarbecovirus monoclonal antibody (mAb) S309 (the parent mAb of sotrovimab) and to the human ACE2 receptor. We provide a blueprint for understanding the marked reduction of binding of other therapeutic mAbs that leads to dampened neutralizing activity. Remodeling of interactions between the Omicron receptor-binding domain and human ACE2 likely explains the enhanced affinity for the host receptor relative to the ancestral virus.

Predicting spike protein NTD mutations of SARS-CoV-2 causing immune evasion by molecular dynamics simulations.

The coronavirus disease 2019 (COVID-19) pandemic was caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Among all the potential targets studied for developing drugs and antibodies, the spike (S) protein is the most striking one, which is on the surface of the virus. In contrast with the intensively investigated immunodominant receptor-binding domain (RBD) of the protein, little is known about the neutralizing antibody binding mechanisms of the N-terminal domain (NTD), let alone the effects of NTD mutations on antibody binding and thereby the risk of immune evasion. Based on 400 ns molecular dynamics simulation for 11 NTD-antibody complexes together with other computational approaches in this study, we investigated critical residues for NTD-antibody binding and their detailed mechanisms. The results show that 36 residues on the NTD including R246, Y144, K147, Y248, L249 and P251 are critically involved in the direct interaction of the NTD with many monoclonal antibodies (mAbs), indicating that the viruses harboring these residue mutations may have a high risk of immune evasion. Binding free energy calculations and an interaction mechanism study reveal that R246I, which is present in the Beta (B.1.351/501Y.V2) variant, may have various impacts on current NTD antibodies through abolishing the hydrogen bonds and electrostatic interaction with the antibodies or affecting other interface residues. Therefore, special attention should be paid to the mutations of these key residues in future antibody and vaccine design and development.

Insights on the mutational landscape of the SARS-CoV-2 Omicron variant receptor-binding domain.

The Omicron variant features enhanced transmissibility and antibody escape. Here, we describe the Omicron receptor-binding domain (RBD) mutational landscape using amino acid interaction (AAI) networks, which are well suited for interrogating constellations of mutations that function in an epistatic manner. Using AAI, we map Omicron mutations directly and indirectly driving increased escape breadth and depth in class 1-4 antibody epitopes. Further, we present epitope networks for authorized therapeutic antibodies and assess perturbations to each antibody's epitope. Since our initial modeling following the identification of Omicron, these predictions have been realized by experimental findings of Omicron neutralization escape from therapeutic antibodies ADG20, AZD8895, and AZD1061. Importantly, the AAI predicted escape resulting from indirect epitope perturbations was not captured by previous sequence or point mutation analyses. Finally, for several Omicron RBD mutations, we find evidence for a plausible role in enhanced transmissibility via disruption of RBD-down conformational stability at the RBDdown-RBDdown interface.

A vaccine built from potential immunogenic pieces derived from the SARS-CoV-2 spike glycoprotein: A computational approximation.

Coronavirus Disease 2019 (COVID-19) represents a new global threat demanding a multidisciplinary effort to fight its etiological agent-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In this regard, immunoinformatics may aid to predict prominent immunogenic regions from critical SARS-CoV-2 structural proteins, such as the spike (S) glycoprotein, for their use in prophylactic or therapeutic interventions against this highly pathogenic betacoronavirus. Accordingly, in this study, an integrated immunoinformatics approach was applied to identify cytotoxic T cell (CTC), T helper cell (THC), and Linear B cell (BC) epitopes from the S glycoprotein in an attempt to design a high-quality multi-epitope vaccine. The best CTC, THC, and BC epitopes showed high viral antigenicity and lack of allergenic or toxic residues, as well as CTC and THC epitopes showed suitable interactions with HLA class I (HLA-I) and HLA class II (HLA-II) molecules, respectively. Remarkably, SARS-CoV-2 receptor-binding domain (RBD) and its receptor-binding motif (RBM) harbour several potential epitopes. The structure prediction, refinement, and validation data indicate that the multi-epitope vaccine has an appropriate conformation and stability. Four conformational epitopes and an efficient binding between Toll-like receptor 4 (TLR4) and the vaccine model were observed. Importantly, the population coverage analysis showed that the multi-epitope vaccine could be used globally. Notably, computer-based simulations suggest that the vaccine model has a robust potential to evoke and maximize both immune effector responses and immunological memory to SARS-CoV-2. Further research is needed to accomplish with the mandatory international guidelines for human vaccine formulations.

Is amyloid fibrillation related to 3D domain swapping for the C-terminal domain of SARS-CoV main protease?

The C-terminal domain of SARS-CoV main protease (Mpro-C) can form 3D domain-swapped dimer by exchanging the α1-helices fully buried inside the protein hydrophobic core, under non-denaturing conditions. Here, we report that Mpro-C can also form amyloid fibrils under the 3D domain-swappable conditions in vitro, and the fibrils are not formed through runaway/propagated domain swapping. It is found that there are positive correlations between the rates of domain swapping dimerization and amyloid fibrillation at different temperatures, and for different mutants. However, some Mpro-C mutants incapable of 3D domain swapping can still form amyloid fibrils, indicating that 3D domain swapping is not essential for amyloid fibrillation. Furthermore, NMR H/D exchange data and molecular dynamics simulation results suggest that the protofibril core region tends to unpack at the early stage of 3D domain swapping, so that the amyloid fibrillation can proceed during the 3D domain swapping process. We propose that 3D domain swapping makes it possible for the unpacking of the amyloidogenic fragment of the protein and thus accelerates the amyloid fibrillation process kinetically, which explains the well-documented correlations between amyloid fibrillation and 3D domain swapping observed in many proteins.

Emergence of novel combinations of SARS-CoV-2 spike receptor binding domain variants in Senegal.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages that carry mutations in the spike gene are of concern for potential impact to treatment and prevention efforts. To monitor for new SARS-CoV-2 mutations, a panel of specimens were sequenced from both wave one (N = 96), and wave two (N = 117) of the pandemic in Senegal by whole genome next generation sequencing. Amongst these genomes, new combinations of SARS-CoV-2 spike mutations were identified, with E484K + N501T, L452R + N501Y, and L452M + S477N exclusively found in second wave specimens. These sequences are evidence of local diversification over the course of the pandemic and parallel evolution of escape mutations in different lineages.

Hotspots for mutations in the SARS-CoV-2 spike glycoprotein: a correspondence analysis.

Spike glycoprotein (Sgp) is liable for binding of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to the host receptors. Since Sgp is the main target for vaccine and drug designing, elucidating its mutation pattern could help in this regard. This study is aimed at investigating the correspondence of specific residues to the SgpSARS-CoV-2 functionality by explorative interpretation of sequence alignments. Centrality analysis of the Sgp dissects the importance of these residues in the interaction network of the RBD-ACE2 (receptor-binding domain) complex and furin cleavage site. Correspondence of RBD to threonine500 and asparagine501 and furin cleavage site to glutamine675, glutamine677, threonine678, and alanine684 was observed; all residues are exactly located at the interaction interfaces. The harmonious location of residues dictates the RBD binding property and the flexibility, hydrophobicity, and accessibility of the furin cleavage site. These species-specific residues can be assumed as real targets of evolution, while other substitutions tend to support them. Moreover, all these residues are parts of experimentally identified epitopes. Therefore, their substitution may affect vaccine efficacy. Higher rate of RBD maintenance than furin cleavage site was predicted. The accumulation of substitutions reinforces the probability of the multi-host circulation of the virus and emphasizes the enduring evolutionary events.

Mutational landscape of SARS-CoV-2 genome in Turkey and impact of mutations on spike protein structure.

The Coronavirus Disease 2019 (COVID-19) was declared a pandemic in March 2020 by the World Health Organization (WHO). As of May 25th, 2021 there were 2.059.941 SARS-COV2 genome sequences that have been submitted to the GISAID database, with numerous variations. Here, we aim to analyze the SARS-CoV-2 genome data submitted to the GISAID database from Turkey and to determine the variant and clade distributions by the end of May 2021, in accordance with their appearance timeline. We compared these findings to USA, Europe, and Asia data as well. We have also evaluated the effects of spike protein variations, detected in a group of genome sequences of 13 patients who applied to our clinic, by using 3D modeling algorithms. For this purpose, we analyzed 4607 SARS-CoV-2 genome sequences submitted by different lab centers from Turkey to the GISAID database between March 2020 and May 2021. Described mutations were also introduced in silico to the spike protein structure to analyze their isolated impacts on the protein structure. The most abundant clade was GR followed by G, GH, and GRY and we did not detect any V clade. The most common variant was B.1, followed by B.1.1, and the UK variant, B.1.1.7. Our results clearly show a concordance between the variant distributions, the number of cases, and the timelines of different variant accumulations in Turkey. The 3D simulations indicate an increase in the surface hydrophilicity of the reference spike protein and the detected mutations. There was less surface hydrophilicity increase in the Asp614Gly mutation, which exhibits a more compact conformation around the ACE-2 receptor binding domain region, rendering the structure in a "down" conformation. Our genomic findings can help to model vaccination programs and protein modeling may lead to different approaches for COVID-19 treatment strategies.

Mechanisms of SARS-CoV-2 Evolution Revealing Vaccine-Resistant Mutations in Europe and America.

The importance of understanding SARS-CoV-2 evolution cannot be overlooked. Recent studies confirm that natural selection is the dominating mechanism of SARS-CoV-2 evolution, which favors mutations that strengthen viral infectivity. Here, we demonstrate that vaccine-breakthrough or antibody-resistant mutations provide a new mechanism of viral evolution. Specifically, vaccine-resistant mutation Y449S in the spike (S) protein receptor-binding domain, which occurred in co-mutations Y449S and N501Y, has reduced infectivity compared to that of the original SARS-CoV-2 but can disrupt existing antibodies that neutralize the virus. By tracking the evolutionary trajectories of vaccine-resistant mutations in more than 2.2 million SARS-CoV-2 genomes, we reveal that the occurrence and frequency of vaccine-resistant mutations correlate strongly with the vaccination rates in Europe and America. We anticipate that as a complementary transmission pathway, vaccine-breakthrough or antibody-resistant mutations, like those in Omicron, will become a dominating mechanism of SARS-CoV-2 evolution when most of the world's population is either vaccinated or infected. Our study sheds light on SARS-CoV-2 evolution and transmission and enables the design of the next-generation mutation-proof vaccines and antibody drugs.

A Routine Sanger Sequencing Target Specific Mutation Assay for SARS-CoV-2 Variants of Concern and Interest.

As SARS-CoV-2 continues to spread among human populations, genetic changes occur and accumulate in the circulating virus. Some of these genetic changes have caused amino acid mutations, including deletions, which may have a potential impact on critical SARS-CoV-2 countermeasures, including vaccines, therapeutics, and diagnostics. Considerable efforts have been made to categorize the amino acid mutations of the angiotensin-converting enzyme 2 (ACE2) receptor binding domain (RBD) of the spike (S) protein, along with certain mutations in other regions within the S protein as specific variants, in an attempt to study the relationship between these mutations and the biological behavior of the virus. However, the currently used whole genome sequencing surveillance technologies can test only a small fraction of the positive specimens with high viral loads and often generate uncertainties in nucleic acid sequencing that needs additional verification for precision determination of mutations. This article introduces a generic protocol to routinely sequence a 437-bp nested RT-PCR cDNA amplicon of the ACE2 RBD and a 490-bp nested RT-PCR cDNA amplicon of the N-terminal domain (NTD) of the S gene for detection of the amino acid mutations needed for accurate determination of all variants of concern and variants of interest according to the definitions published by the U.S. Centers for Disease Control and Prevention. This protocol was able to amplify both nucleic acid targets into cDNA amplicons to be used as templates for Sanger sequencing on all 16 clinical specimens that were positive for SARS-CoV-2.

Immunizations with diverse sarbecovirus receptor-binding domains elicit SARS-CoV-2 neutralizing antibodies against a conserved site of vulnerability.

Viral mutations are an emerging concern in reducing SARS-CoV-2 vaccination efficacy. Second-generation vaccines will need to elicit neutralizing antibodies against sites that are evolutionarily conserved across the sarbecovirus subgenus. Here, we immunized mice containing a human antibody repertoire with diverse sarbecovirus receptor-binding domains (RBDs) to identify antibodies targeting conserved sites of vulnerability. Antibodies with broad reactivity against diverse clade B RBDs targeting the conserved class 4 epitope, with recurring IGHV/IGKV pairs, were readily elicited but were non-neutralizing. However, rare class 4 antibodies binding this conserved RBD supersite showed potent neutralization of SARS-CoV-2 and all variants of concern. Structural analysis revealed that the neutralizing ability of cross-reactive antibodies was reserved only for those with an elongated CDRH3 that extends the antiparallel beta-sheet RBD core and orients the antibody light chain to obstruct ACE2-RBD interactions. These results identify a structurally defined pathway for vaccine strategies eliciting escape-resistant SARS-CoV-2 neutralizing antibodies.

SARS-CoV-2 Variants, RBD Mutations, Binding Affinity, and Antibody Escape.

Since 2020, the receptor-binding domain (RBD) of the spike protein of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been constantly mutating, producing most of the notable missense mutations in the context of "variants of concern", probably in response to the vaccine-driven alteration of immune profiles of the human population. The Delta variant, in particular, has become the most prevalent variant of the epidemic, and it is spreading in countries with the highest vaccination rates, causing the world to face the risk of a new wave of the contagion. Understanding the physical mechanism responsible for the mutation-induced changes in the RBD's binding affinity, its transmissibility, and its capacity to escape vaccine-induced immunity is the "urgent challenge" in the development of preventive measures, vaccines, and therapeutic antibodies against the coronavirus disease 2019 (COVID-19) pandemic. In this study, entropy-enthalpy compensation and the Gibbs free energy change were used to analyze the impact of the RBD mutations on the binding affinity of SARS-CoV-2 variants with the receptor angiotensin converting enzyme 2 (ACE2) and existing antibodies. Through the analysis, we found that the existing mutations have already covered almost all possible detrimental mutations that could result in an increase of transmissibility, and that a possible mutation in amino-acid position 498 of the RBD can potentially enhance its binding affinity. A new calculation method for the binding energies of protein-protein complexes is proposed based on the entropy-enthalpy compensation rule. All known structures of RBD-antibody complexes and the RBD-ACE2 complex comply with the entropy-enthalpy compensation rule in providing the driving force behind the spontaneous protein-protein docking. The variant-induced risk of breakthrough infections in vaccinated people is attributed to the L452R mutation's reduction of the binding affinity of many antibodies. Mutations reversing the hydrophobic or hydrophilic performance of residues in the spike RBD potentially cause breakthrough infections of coronaviruses due to the changes in geometric complementarity in the entropy-enthalpy compensations between antibodies and the virus at the binding sites.

Specific delivering of RNAi using Spike's aptamer-functionalized lipid nanoparticles for targeting SARS-CoV-2: A strong anti-Covid drug in a clinical case study.

Coronavirus (SARS-CoV-2) as a global pandemic has attracted the attention of many scientific centers to find the right treatment. We expressed and purified the recombinant receptor-binding domain (RBD) of the SARS-CoV-2 spike (S) protein, and specific RBD aptamers were designed using SELEX method. RNAi targeting nucleocapsid phosphoprotein was synthesized and human lung cells were inoculated with aptamer-functionalized lipid nanoparticles (LNPs) containing RNAi. The results demonstrated that RBD aptamer having KD values of 0.290 nm possessed good affinity. Based on molecular docking and efficacy prediction analysis, siRNA molecule was showed the best action. LNPs were appropriately functionalized by aptamer and contained RNAi molecules. Antiviral assay using q-PCR and ELISA demonstrated that LNP functionalized with 35 µm Apt and containing 30 nm RNAi/ml of cell culture had the best antiviral activity compared to other concentrations. Applied aptamer in the nanocarrier has two important functions. First, it can deliver the drug (RNAi) to the surface of epithelial cells. Second, by binding to the SARS-CoV-2 spike protein, it inhibits the virus entrance into cells. Our data reveal an interaction between the aptamer and the virus, and RNAi targeted the virus RNA. CT scan and the clinical laboratory tests in a clinical case study, a 36-year old man who presented with severe SARS-CoV-2, demonstrated that inhalation of 10 mg Apt-LNPs-RNAi nebulized/day for six days resulted in an improvement in consolidation and ground-glass opacity in lungs on the sixth day of treatment. Our findings suggest the treatment of SARS-CoV-2 infection through inhalation of Aptamer-LNPs-RNAi.

α-Hemolysin-Aided Oligomerization of the Spike Protein RBD Resulted in Improved Immunogenicity and Neutralization Against SARS-CoV-2 Variants.

The increase in confirmed COVID-19 cases and SARS-CoV-2 variants calls for the development of safe and broad cross-protective vaccines. The RBD of the spike protein was considered to be a safe and effective candidate antigen. However, the low immunogenicity limited its application in vaccine development. Herein, we designed and obtained an RBD heptamer (mHla-RBD) based on a carrier protein-aided assembly strategy. The molecular weight of mHla-RBD is up to 450 kDa, approximately 10 times higher than that of the RBD monomer. When formulated with alum adjuvant, mHla-RBD immunization significantly increased the immunogenicity of RBD, as indicated by increased titers of RBD-specific antibodies, neutralizing antibodies, Th2 cellular immune response, and pseudovirus neutralization activity, when compared to RBD monomer. Furthermore, we confirmed that RBD-specific antibodies predominantly target conformational epitopes, which was approximately 200 times that targeting linear epitopes. Finally, a pseudovirus neutralization assay revealed that neutralizing antibodies induced by mHla-RBD against different SARS-CoV-2 variants were comparable to those against the wild-type virus and showed broad-spectrum neutralizing activity toward different SARS-CoV-2 variants. Our results demonstrated that mHla-RBD is a promising candidate antigen for development of SARS-CoV-2 vaccines and the mHla could serve as a universal carrier protein for antigen design.